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ABSTRACT 

The discovery of the Cosmic Infrared Background (GIB) in 1996, together with 
recent cosmological surveys from the mid-infrared to the miUimeter have revolu- 
tionized our view of star formation at high redshifts. It has become clear, in the 
last decade, that a population of galaxies that radiate most of their power in the 
far- infrared (the so-called "infrared galaxies" ) contributes an important part of the 
whole galaxy build-up in the Universe. Since 1996, detailed (and often painful) 
investigations of the high-redshift infrared galaxies have resulted in the spectacular 
progress covered in this review. We outline the nature of the sources of the GIB 
including their star-formation rate, stellar and total mass, morphology, metallicity 
and clustering properties. We discuss their contribution to the stellar content of 
the Universe and their origin in the framework of the hierarchical growth of struc- 
tures. We finally discuss open questions for a scenario of their evolution up to the 
present-day galaxies. 

Subject headings: cosmology - evolution - luminosity function - starburst - star 
formation 



1. INTRODUCTION 

The Cosmic Infrared Background (GIB) can 
be defined as the part of the present radia- 
tion content of the Universe that is made es- 
sentially of the long wavelength output from all 
sources throughout the history of the Universe. 
The radiation content in the microwave part of 
the spectrum is dominated by the Cosmic Mi- 
crowave Background (CMB) produced in the 
hot and dense phases of the universe. It domi- 
nates for frequencies below 800 GHz. Neverthe- 
less the very different spectra of the GIB with 
respect to the CMB (both its purely Planckian 
part and its Compton distortion expected to be 
the dominant one at these frequencies) allow 
them to be separated very efhciently down to 
frequencies close to the peak of the CMB (150 
GHz). Furthermore in this frequency regime 
the GIB dominates the galactic emission in the 
lowest cirrus regions by a factor ~4. The cos- 
mic background due to sources (CMB excluded) 
presents two maxima: one in the optical, one in 
the far-infrared, with roughly equal brightness 
(in vli,) and with a minimum around 5 /zm. 
This minimum is created by the decrease of 
brightness of the stellar component with wave- 
length combined with the rising brightness of 



the dust, very small grains, and of the Active 
Galactic Nuclei (AGN) non thermal emission in 
the thermal infrared. The GIB is defined as the 
cosmic background at wavelengths longward of 
this minimum. An understanding of the na- 
ture and redshift distribution of the sources of 
the GIB, although relatively new, is an integral 
part of the understanding of the formation and 
evolution of galaxies. 

The standard hierarchical model of structure 
formation has received strong observational 
support from observations of the large-scale 
distribution of galaxies, clusters, intergalac- 
tic clouds, combined with CMB anisotropies 
that constrain the initial large scale power 
spectrum within the concordance cosmological 
model framework. Scenarios for galaxy forma- 
tion and evolution can be confronted with the 
very quickly rising set of observations of ex- 
tragalactic sources at higher and higher red- 
shifts. Nevertheless many critical questions re- 
main open on the cooling of collapsed struc- 
tures, angular momentum of galaxies, star for- 
mation efficiency. Initial Mass Function (IMF) 
of the stars formed, role of feedback mecha- 
nisms, the physics and role of merging and ac- 
cretion in the construction of galaxies. 
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As ultraluminous infrared galaxies were found 
to be often associated with mergers or interact- 
ing galaxies, it can be expected that the sources 
of the CIB carry critical information about the 
history of merging (e.g., Sanders & Mirabel 
1996: Genzcl & Ccsarsky 2000). Because about 
half of the energy from extragalactic sources is 
in the CIB, the determination of the source of 
this energy (starbursts or massive black hole ac- 
cretion in dust enshrouded AGNs) should shed 
some light on how galaxies evolve. 

Since the discovery of the CIB in the COBE 
data (e.g., Hauser & Dwek 2001), identify- 
ing the sources of the CIB, their redshift dis- 
tribution, and nature progressed at increasing 
speed especially through multiwavelength anal- 
ysis. This review attempts to give the broad- 
band observational picture for the identification 
of the sources of the CIB. The detailed analy- 
sis of individual infrared galaxies is outside the 
scope of this review. We then discuss implica- 
tions, both well-established ones and tentative 
ones, as well as directions for future work. The 
spectroscopy aspects are covered in the review 
by Solomon & Vanden Bout (2005, in this vol- 
ume) . This is at present a fascinating but mov- 
ing target! The coming decade will be a very 
rich one. New, very powerful long-wavelengths 
observation tools will bring many striking new 
results like the Spitzer observatory or the be- 
ing built experiments like Herschel and ALMA. 
Throughout this review, the cosmology is fixed 
to Oa=0.7, ^lm=0.3, Ho=100 h km s'^ Mpc'^ 
with h=0.65. 

2. DUST IN THE LOCAL UNIVERSE 

A fraction of the stellar radiation produced 
in galaxies is absorbed by dust and re-radiated 
from mid-infrared to millimeter wavelengths. 
Understanding dust properties and the associ- 
ated physics of the absorption and emission are 
thus essential. These determine the Spectral 
Energy Distribution (SED) of the galaxies. 

2.1. Dust Particles 

Small dust particles with sizes ranging from 
a nanometer to a fraction of micrometer are 
ubiquitous in the interstellar medium. They 
result from natural condensation in cool stel- 
lar atmospheres, supernovae, and the interstel- 
lar medium of the heavy elements produced 
by the nucleosynthesis in stars and released to 



the diffuse medium by late type stars and su- 
pernovae explosions. Interstellar grain models 
have been improved for 30 years in order to fit 
all observational constraints: elemental abun- 
dances of the heavy elements, UV, visible and 
infrared absorption and scattering properties, 
infrared emission, polarization properties of the 
absorbed and emitted light. The models in- 
clude a mixture of amorphous silicate grains 
and carbonaceous grains, each with a wide size 
distribution ranging from molecules containing 
tens of atoms to large grains > 0.1 /im in di- 
ameter that can be coated with ices in dense 
clouds and/or organic residues (e.g.. Desert et 
al. 1990; Li & Drainc 2001). It is now widely 
accepted that the smallest carbonaceous grains 
are Polycyclic Aromatic Hydrocarbons (PAHs) 
that emit a substantial fraction of the energy 
in a set of features between 3 and 17 /im (3.3, 
6.2, 7.7, 8.6, 11.3, 12.7, 16.3, 17 /xm for the 
main ones) that used to be known as the UIB 
for Unidentified Infrared Bands. These features 
result from C-C and C-H stretching/bending 
vibrational bands excited by the absorption of 
a single UV or optical photon and are a good 
tracer of normal and moderately active star for- 
mation activity in spiral and irregular galaxies 
(e.g., Helou et al. 2000; Peeters et al. 2004). 
For radii a>50 A, the carbonaceous grains are 
often assumed to have graphitic properties. The 
so-called very small grains of the interstellar 
medium are small enough to have very low heat 
capacity, so their temperature are significantly 
affected by single-photon absorption. In the 
diffuse ISM of our Galaxy, they dominate the 
infrared emission for wavelengths smaller than 
about 80 /im. At longer wavelengths, the in- 
frared spectrum is dominated by the emission 
of the larger grains at their equilibrium temper- 
atures. Considering the energy density of the 
radiation in a galactic disc like ours, the tem- 
perature of the larger grains is rather low; 15 to 
25 K. For these grains, the far- infrared emissiv- 
ity decreases roughly as the square of the wave- 
length. This in turn makes the temperature 
dependence on the radiation energy density u 
very weak (T ~ u^^^). For a galaxy like the 
Milky Way, the infrared part of the SED peaks 
at 170 /um whereas for an Ultra Luminous In- 
frared Galaxy (ULIRG) it peaks at about 60 
/xm: a factor 3 in temperature for a factor 10'^ 
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Fig. 1. — Spectral energy distributions of galaxies from UV to the millimeter. The ULIRG is observed at redshift 
;=0.66 and is represented here in the rest-frame (from Galliano 2004). 



in energy density or luminosity. At long wave- 
lengths in the submillimeter and millimeter, the 
intensity should increase like 7,^ ~ i^^. 



2.2. Extinction 

In our Galaxy, the extinction curve of the dif- 
fuse ISM has been known for a few decades. The 
average optical depth perpendicular to the disk 
of our Galaxy in the solar vicinity is small {Ay 
~ 0.2) and typical of spiral galaxies. The av- 
erage optical depth increases to a few in large 
molecular cloud complexes. It can become very 
large in galactic nuclei. Finally it should also 
be remembered that the optical depth in the 
UV is typically 10 times larger than that in op- 
tical wavelengths. The conversion of star light 
into infrared radiation will thus depend strongly 
on the location of the stars and their spectral 
types. 

In external galaxies, modeling the extinction 



is very hard because it strongly depends on 
the geometric distribution of the ISM and of 
the chemical abundances. Simple models have 
been used to take this into account to first or- 
der. Galaxies can be modeled as an oblate ellip- 
soid where absorbers (dust) and sources (stars) 
are homogeneously mixed; the dust absorption 
can be computed in a "screen" or "sandwich" 
geometry (dust layers in front of the stars or 
sandwiched between two star layers). As a con- 
sequence, the reddening curve average over a 
whole galaxy appears to vary within a class of 
objects and between the different classes, from 
normal star-forming galaxies to highly concen- 
trated starburst. It is thus very difficult to de- 
rive the total dust optical depth (e.g., Calzetti 
et al. 1994). In the local Universe, the average 
extinction per galaxy is quite low. About one 
third of the bolometric flux is emitted in the 
far-infrared, and this is typical of our Galaxy. 
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In more actively star-forming galaxies, up to 
70% of the bolometric flux is emitted in the far- 
infrared. In some of these, the starburst activity 
is mostly in the disk (like in M51). For a given 
total luminosity, the radiation energy density is 
lower than in the case of a starburst concen- 
trated in a small volume in the nucleus. In this 
case, the dust will be hotter due to the larger 
energy density and the conversion of stellar light 
to infrared will be more efficient. Some ULIRGs 
emit more than 95% of their energy in the far- 
infrared (e.g., Arp 220). Such galaxies are very 
compact, dusty starbursts where dust optical 
depths are very large. In such galaxies, fine 
structure and recombination line ratios imply 
an equivalent "screen" dust extinction between 
Av~5 and 50. The result is that the SED is 
significantly distorted in the opposite way from 
the higher dust temperature (less mid-infrared 
emission). In the following, we will refer to "in- 
frared galaxies" and to "optical galaxies" to des- 
ignate galaxies in which the infrared emission, 
respectively optical emission dominates. Differ- 
ent typical spectra of galaxies are shown in Fig- 
ure^from the UV to the millimeter. We clearly 
see the variation of the optical to infrared en- 
ergy ratio as starburst activity increases. 

2.3. Local Infrared Galaxies 

A few very luminous infrared galaxies were 
observed in the seventies (Rieke & Lebof- 
sky 1979). Then IRAS sateUite, launched in 
1983 gave for the first time a proper cen- 
sus of the infrared emission of galaxies at low 
redshift. The Luminosity Function (LF) at 
60 and 100 /im is dominated by Li, spiral 
galaxies as could be expected - the reradi- 
ated stellar luminosity absorbed by dust. In 
addition, a high-luminosity tail of luminous 
galaxies was found (e.g., Sanders & Mirabel 
1996). This high- luminosity tail can be ap- 
proximated by a power-law, ^{L) cx L^^^, 
which gives a space density for the most 
luminous infrared sources well in excess of 
predictions based on the optical LF. These 
sources comprise the Luminous Infrared Galax- 
ies, LIRGs, and the ULIRGs with luminosities 
ll<log(LiR/Lo)<12 and log(LiR,/L0)>12, re- 
spectively. These galaxies are often associated 
with interacting or merging, gas-rich disks. The 
fraction of strongly interacting/merger systems 
increases from ~10% at log(LiR/L0)=lO.5-ll 



to ~100% at log(LiR/Lo)>12. LIRGs are the 
site of intense starburst activity (about 10- 
100 Mq year~^) induced by the interaction 
and/or strong spiral structure. The ULIRG 
phase occurs near the end of the merging pro- 
cess when the two disks overlap. Such galax- 
ies may be the precursors of Quasi Stellar Ob- 
jects (QSOs; Sanders et al. 1988a, 1988b; 
Veilleux et al. 1995; Lutz et al. 1999). 
These objects have been the subject of in- 
tense debate concerning the nature of the dom- 
inant source of emission: starburst versus dust- 
enshrouded AGN (e.g., Filipenko 1992; Sanders 
& Mirabel 1996; Joseph 1999). Indeed, spec- 
tra show evidence of extremely large optical 
depth (heavily reddened continuum and large 
Balmer decrement) but also exhibit AGN-like 
high excitation fine-structure lines. We had to 
wait for ISO to clearly determine the power 
sources of ULIRGs. The difference between 
the mid-infrared spectra of starburst and AGNs 
is striking. Starburst are often characterized 
by strong, low-excitation, fine-structure lines, 
prominent PAH features and a weak A >10 /xm 
continuum whereas AGNs display a highly ex- 
cited emission line spectrum with weak or no 
PAH features, plus a strong mid-infrared con- 
tinuum. It has been thus possible to build mid- 
infrared diagnostic diagrams (e.g., Genzel et al. 
1998; Laurent et al. 2000) that clearly sepa- 
rates starburst-dominated galaxies from AGN- 
dominated galaxies. These diagrams demon- 
strate that ULIRGs appear to be composite ob- 
jects, but star formation dominates in most ob- 
jects. That is on average, >70% of the rera- 
diated energy comes from starbursts and < 
30% comes from AGNs (Genzel et al. 1998; 
Lutz et al. 1998). However the fraction of 
AGN-powered objects increases with luminos- 
ity. About 15% of ULIRGs at luminosities be- 
low 2x10^^ Lq are mostly AGN powered, but 
this fraction increases to about half at higher 
luminosity. 

All these well-studied LIRGs and ULIRGs are 
at low redshift. They do not dominate the en- 
ergy production locally. As an example, the 
total infrared luminosity from these galaxies in 
the IRAS Bright Galaxy Sample accounts for 
only ^6% of the infrared emission in the local 
Universe (Soifer & Neugeubauer, 1991). As we 
will see, the situation changes dramatically at 
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higher redshift where these galaxies fully dom- 
inate the infrared energy output. 

3. THE COSMIC INFRARED BACKGROUND 

The CIB is the infrared part of the extra- 
galactic background, the radiation content of 
the Universe today produced by galaxies at all 
redshifts and seen as an isotropic extragalac- 
tic background radiation. Patridge & Peebles 
(1967) predicted that observations of such a 
background could give powerful constraints on 
the cosmological evolution. 

3.1. General Observations and Direct 
Cosmological Implications 

The detection of the infrared part of the 
extragalactic background (the CIB for Cosmic 
Infrared Background) was the major objective 
of the DIRBE experiment aboard COBE. 
In fact, the CIB was first detected at long 
wavelengths by using the FIRAS spectrometer: 
A > 200 /im (Puget et al. 1996). The CIB 
has subsequently been detected by DIRBE 
at 2.4, 3.5, 100, 140, 240 ^m (see Hauser 
& Dwek 2001 and Kashlinsky 2005 for two 
reviews). The extragalactic background at 2.4 
and 3.5 /im is significantly larger than that 
predicted by the integrated galaxy counts and 
their extrapolation. Similarly, the extragalactic 
background in the optical has been finally 
evaluated in combining several methods by 
Bernstein ct al. (2002) and found to be larger 
than the value given by the integrated fluxes 
of galaxies by a factor larger than 2. In the 
mid-infrared, the interplanetary zodiacal dust 
emission is so strong that only upper limits 
were obtained by DIRBE. The combination 
of number counts by ISO/ISOCAM at 15 /im 
(see Elbaz & Cesarsky 2003) and by Spitzer at 
24 /xm (e.g., Papovich et al. 2004) giving lower 
limits, with the observations of TeV gamma 
ray emission from distant AGNs (e.g., Renault 
et al. 2001; Dwek & Krennrich 2005), gives 
a good measurement of the background at 
these wavelengths. The full cosmic background 
spectrum is shown in Figure |21 Only most 
recent and strongly constraining measurements 
have been plotted for clarity. 

Figure |21 clearly shows that the optical and 
infrared cosmic backgrounds are well separated. 
The first surprising result is that the power in 



the infrared is comparable to the power in the 
optical. In contrast, we know that locally, the 
infrared output of galaxies is only one third of 
the optical output. This implies that infrared 
galaxies grow more luminous with increasing z 
faster than do optical galaxies. A second im- 
portant property to note is that the slope of 
the long wavelength part of the CIB, oc z^^-* 
(Gispert et al. 2000), is much less steep than 
the long wavelengths spectrum of galaxies (as il- 
lustrated in Figure El with the M82 SED). This 
implies that the millimeter CIB is not due to the 
millimeter emission of the galaxies that account 
for the peak of the CIB (~ 150/im). The im- 
plications in terms of energy output have been 
drawn by, e.g. Gispert et al. (2000). The in- 
frared production rate per comoving unit vol- 
ume (a) evolves faster between redshift zero 
and 1 than the optical one and (6) has to stay 
roughly constant at higher redshifts up to red- 
shift 3 at least. 

3.2. The Status of Deep Surveys: Resolved 
Fraction of the >10 fim CIB 

Many surveys from the mid-infrared to the 
millimeter have aimed to resolve the CIB into 
discrete sources. From short to long wave- 
lengths the significant surveys are the following: 

• ISOCAM 15 fim: Three kinds of surveys 
have been done. The shallowest is the 
ELAIS survey (European Large- Area ISO 
Survey, Oliver et al. 2000). The deep- 
est is the survey in the HDF-N (Aussel 
et al. 1999) as well as the surveys in the 
direction of the galaxy clusters (e.g. Met- 
calfe et al. 2003). Altogether, about 1000 
galaxies were detected. Above a sensitiv- 
ity limit of 50 fJ.3y, they produce a 15 /xm 
extragalactic background light of (2.4 ± 
0.5) nW m-2 sr"! (Elbaz et al. 2002). 
This accounts for about 80% of the CIB 
at 15 /im based on the simplest extrapo- 
lation of the counts. 

• Spitzer 24 /im: Spitzer surveys are ongo- 
ing, but more than 10'^ sources have al- 
ready been detected at 24 /im. Integrat- 
ing the first number counts down to 60 
/iJy, a lower limit to the CIB at 24 /tm 
of (1.9±0.6) nW m^^ sr^^ is derived (Pa- 
povich et al. 2004). This accounts for 
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Fig. 2. — The extragalactic background over three decades in frequency from the near UV to milhmeter wavelengths. 
Only strongly constraining measurements have been reported. We show for comparison in grey an SED of M82 
(Chanial, 2003)- a starburst galaxy at L=3x 10^" Lq - normalized to the peak of the CIB at 140 /im. References 
for data points are given in Table [Tl 



TABLE 1 

Extragalactic background references for Figure |^ 



Wavelength (fJ-m) 


Experiment 


Measurement 


Reference 


0.2 


FOCA 


Number counts & model 


Armand et al. 1994 


0.30, 0.56, 0.81 


HST/Las 


Diffuse emission 


Bernstein et al. 2002 




Campanas Obs. 




Mattila 2003 


2.2< A <4 


IRTS 


Diffuse emission 


Matsumoto et al. 2005 


2.2, 3.3 


DIRBE/Lick 


Diffuse emission 


Gorjian et al. 2000 


1.25, 2.2 


DIRBE/2MASS 


Diffuse emission 


Wright 2001 








Cambresy et al. 2001 


10, 15 


CAT 


7-rays 


Renault et al. 2001 


15 


ISO/ISOCAM 


Number counts 


Elbaz et al. 1999 


24 


Spitzer/MlPS 


Number counts 


Papovich et al. 2004 


60 


IRAS 


Power spectrum 


Miville-Deschenes et 








al. 2001 


100 


DIRBE 


Diffuse emission 


Renault et al. 2001 


140, 240 


DIRBE/WHAM 


Diffuse emission 


Lagache et al. 2000 


140, 240 


DIRBE 


Diffuse emission 


Hauser et al. 1998 


850 


SCUBA 


Number counts 


Small et al. 2002 


200< A <1200 


EZRAS 


Diffuse emission 


Lagache et al. 2000 



about 70% of the CIB at 24 iim, based on 



a simple extrapolation of the counts. 
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Fig. 3. — Sensitivity to the bolometric luminosity and star-formation rate, assuming star forming galaxies of various 
infrared and submillimeter experiments. Detections of at least 10 sources in the surveys can be expected in the areas 
above the curves. We assumed the scenario of a typical deep survey (when available). IRAS 60 fim (Sjy > 1 Jy, all sky); 
ISOCAM 15 fim {S^ > 250/iJy, 2 Sq. Deg.); ISOPHOT 170 fj.m (S^ > 180 mjy, 5 Sq. Deg.); Spitzer/MIPS 24 fim 
(S^ > 80/iJy, 5 Sq. Deg.); Spitzer/MIPS 70 ^im (S^ > 25 mJy, 5 Sq. Deg.); Spitzer/MIPS 160 fim (S^ > 50 mJy, 5 
Sq. Deg.); SCUBA 850 fim {Sv > 1 mJy, 1 Sq. Deg.). This plot makes use of the Lagache et al. (2004) model (see 
the Appendix). 



• Spitzer 70 /im: Counts have been derived 
down to 15 mJy. Integrating these counts 
corresponds 0.95 nW sr~^ which ex- 
plains ~20% of the CIB at 70 fim as 
derived from the Lagache et al. (2004) 
model (Dole et al. 2004a). 

• ISOPHOT 90 /im: The most relevant 
data comes from the ELAIS survey that 
covers about 12 square degrees at 90 /xm. 
Counts have been obtained down to 95 
mJy (Heraudeau et al. 2004) resolving 
less than 5% of the CIB as derived from 
DIRBE by Renauh et al. (2001). 

• Spitzer 160 /im: First counts are derived 
down to 50 mJy. The integral of these 
counts corresponds to 1.4 nW m^^ sr~^, 



which explain about 7% of the CIB at 160 
/im (Dole et al. 2004a). 

• ISOPHOT 170 /im: Two main surveys 
have been conducted: the FIRBACK (La- 
gache & Dole 2001) and the Lockman hole 
(Kawara et al. 2004) surveys that cover- 
ing about 5 square degrees. Counts down 
to 135 mJy contribute to less than 5% of 
the CIB (Dole et al. 2001). 

• SCUBA 450 fim: Deep surveys at 450 /im 
are very hard to conduct from the ground. 
A few galaxies are detected between 10 
and 50 mJy in deep surveys. Number 
counts down to 10 mJy give a lower limit 
on the CIB of about (0.7 ± 0.4) nW m'^ 
sr-i (Small et al. 2002). This resolves 
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about 15% of the CIB at 450 fim. 

• SCUBA 850 iim: Over 500 arcmin^ of 
blank sky has been surveyed by several 
groups using SCUBA. The observations 
range from an extremely deep survey in 
the area of the HDF-N to a wider-field 
shallower survey. Also about 40 arcmin^ 
of lensed cluster fields have been observed 
in which 17 sources have been detected. 
Fluxes range from about 0.5 to 8 mJy. 
The flux density in the resolved submil- 
limeter population down to 1 mJy is (0.3 
± 0.1) nW m-2 sr-i (Small et al. 2002). 
This account for 60% of the CIB at 850 
fim. Note that in deep surveys, sources 
with S85o> 3 mJy contribute to '^30% of 
the CIB. 

• MAMBO 1200 fi-ai: Deep blank sky ar- 
eas surveyed by MAMBO cover about 500 
arcmin^. Number counts are given by 
Greve et al. (2004). By integrating the 
counts from 2.25 to 5.75 mJy, the resolved 
CIB is about 0.02 nW m~^ sr~^, or about 
10% of the total CIB at 1200 /im. 

FigureOshows the capabilities of the different 
surveys to find distant LIRGs. Spitzer observa- 
tions at 24 //m are the most powerful tool to find 
LIRGs up to z '-2.2; ISOCAM was limited at 
z ~1.2. Distant ULIRGs are found by deep and 
large surveys at 24 and 850 /im. Note that ca- 
pabilities have been computed using the model 
of Lagache et al. (2004). This empirical model 
is based on only two populations of galaxies; it 
aims only to model the redshift evolution of the 
average population. It reproduces all the obser- 
vations from mid-infrared to the millimeter (see 
Appendix). Lewis et al. (2005) showed that 
a more sophisticated, bivariate SED does not 
much change the average properties although it 
does significantly change the dispersion. The 
Lagache et al. (2004) model is thus used in this 
paper as a tool to discuss observations and pre- 
dictions. 

3.3. Redshift Contribution to the CIB 

From Figure |21 we see that contributions 
from galaxies at various redshifts are needed 
to fill the CIB SED shape. The bulk of the 
CIB in energy, i.e., the peak at about 150 fim, 
is not resolved in individual sources but one 



TABLE 2 

Redshift at which the Cosmic Infrared 
Background is resolved at 20, 50 or 80% (from 
THE model of Lagache et al. 2004). 



Wavelength 


20% 


50% 


80% 


15 /im 


0.5 


1.0 


1.3 


24 /im 


0.5 


1.3 


2.0 


70 /im 


0.5 


1.0 


1.5 


100 fim 


0.7 


1.0 


1.7 


160 iim 


0.7 


1.3 


2.0 


350 iim 


1.0 


2.0 


3.0 


850 


2.0 


3.0 


4.0 


1.4 mm 


2.5 


3.5 


4.5 


2.1 mm 


2.0 


3.5 


5.0 



dominant contribution at the CIB peak can 
be inferred from the ISOCAM deep surveys. 
ISOCAM galaxies with a median redshift 
of ~0.7 resolve about 80% of the CIB at 
15 /im. Elbaz et al. (2002) separate the 15 /zm 
galaxies into different classes (ULIRGs, LIRGs, 
Starbursts, normal galaxies and AGNs) and 
extrapolate the 15 fluxes to 140 /im using 
template SEDs. A total brightness of (16±5) 
nW m~^ sr~^ is found, which makes up about 
two thirds of the CIB observed at 140 /im by 
COBE/DIRBE. Hence, the galaxies detected 
by ISOCAM are responsible for a large fraction 
of energy of the CIB. About one half of the 
140 /im CIB is due to LIRGs and about one 
third to ULIRGs. However, these ISOCAM 
galaxies make little contribution to the CIB 
in the millimeter and submillimeter. There, 
the CIB must be dominated by galaxies at 
rather high redshift for which the SED peak 
has been shifted. The redshift contribution 
to the CIB is illustrated in Figure 0] We 
clearly see that the submillimeter/millimeter 
CIB contains information on the total energy 
output by the high- redshift galaxies {z >2). 
This is supported by the redshift distribution 
of the SCUBA sources at 850 /xm with S85o>3 
mJy that make about 30% of the CIB and have 
a median redshift of 2.2 (Chapman et al. 2005). 

FigureElshows the fraction of resolved CIB as 
a function of redshift for selected wavelengths. 
Fifty percents of the CIB is due to galaxies at 
redshift below 1 at 15 and 70 /im, 1.3 at 24 
and 160 /im, 2 at 350 /im, 3 at 850 /xm and 
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Fig. 4. — Cumulative contribution to the CIB of galaxies at various redshifts from 0.5 to 8, from the model of 
Lagache et al. (2004). Measurements of the CIB are reported with the same symbols as in Figure 1^ 



3.5 at 2000 /im (see also Table I^J. It is clear 
that from the far-infrared to the millimeter, 
the CIB at longer wavelengths probes sources 
at higher redshifts. 

From Section IT^ we see that the most con- 
straining surveys in term of resolving the CIB 
are those at 15, 24 and 850 /zm. Moreover, the 
capabilities of these surveys to find high- 2; ob- 
jects are the best among all other existing sur- 
veys (see Figure 13). These surveys probe the 
CIB in well-defined and distinct redshift ranges, 
with median redshifts of 0.7 (Liang et al. 2004), 
~1 (Caputi et al. 2005 and L. Yan, private com- 
munication), and 2.2 (Chapman et al. 2005) at 
15, 24 and 850 /zm, respectively. Such well- 
defined redshift ranges can be understood by 
looking at the K-correction. The K-correction 



is defined as: 



KiL,z) 



(1) 



L,{z = 0) 

where iy(z = 0) is the rest- frame luminosity. 
This correction is specific to the spectrum of 
the population considered at a given luminosity 
and redshift. FigureElshows the K-correction at 
15, 24, 70, 160 and 850 fim. The broad plateau 
observed around z = 1 at 15 fim and around 
z = 2 at 24 /im is caused by the PAHs' fea- 
tures. At longer wavelengths, the slow decrease 
of the K-correction is caused by the shape of 
the starburst spectra around the peak of their 
emission. At 850 /im, the monotonic rise favors 
the detection of high-z objects. 

4. GALAXIES AT REDSHIFTS 0.5 < Z < 1.5 

At the time this review was written, most of 
the detailed informations on dusty galaxies in 
the 0.5 < z < 1.5 redshift range comes from 
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Fig. 5. — Cumulative fraction of the CIB content as a function of redshift for various wavelengths, from the model 
of Lagache et al. (2004). 



galaxies selected with the ISOCAM cosmologi- 
cal surveys and the multi- wavelength analysis of 
detected sources. The ISOCAM extragalactic 
surveys were performed with two filters, LW2 
(5-8.5 /im) and LW3 (12-18 /xm) centered at 
6.75 and 15 /xm, respectively. However, be- 
cause of star contamination and because the 
stellar light dominates the flux from galaxies 
in the 6.75 fim band above redshift 0.4, only 
the 15 fim surveys are relevant here. Shallow, 
deep and ultra-deep surveys were performed 
in various fields including the Lockman hole, 
Marano, northern and southern Hubble Deep 
Field (HDF), and Canada- France Redshift Sur- 
vey (CFRS) (e.g., Aussel et al. 1999; Flores et 
al. 1999; Lariet al. 2001; Gruppioni et al. 2002; 
Mann ct al. 2002; Elbaz & Cesarsky 2003; Sato 
et al. 2003). Deeper images have been made in 
the direction of distant clusters (e.g., Metcalfe 
et al. 2003). Finally, the bright end of the lu- 



minosity function was explored by the ELAIS 
survey (e.g., Oliver et al. 2000). The deep- 
est surveys reach a completeness limit of about 
100 /LtJy at 15 /zm (without lensing). The most 
relevant data to this section are the deep and 
ultra-deep surveys. 

4.1. Detailed Properties 

To find out the nature and redshift distribu- 
tion of the 15 /im deep survey sources, many 
foUowup observations have been conducted in- 
cluding HST imaging and VLT spectroscopy. 
With a point-spread function full width at half 
of maximum of 4.6 arcsec at 15 ^m, optical 
counterparts are easily identified. Redshifts are 
found using emission and/or absorption lines. 
From field to field, the median redshift varies 
from 0.52 to 0.8, a quite large variation due 
to sample variance. Each field clearly exhibits 
one or two redshift peaks, with velocity disper- 
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Fig. 6. — K-correction at 15, 24, 70 and 160 and 850 f-tm for a typical LIRG with L=2xl0" Lq (fr om the model 
of Lagache ot al. 2004). 



sion characteristic of clusters or galaxy groups. 
Most of ISO CAM galaxies have redshifts be- 
tween '^0.3 and 1.2, consistent with Figure |3| 
About 85% of the ISO galaxies show obvious 
strong emission hues (e.g., [Oil] 3723, H.^, H^, 
[OIII] 4959, 5007). These hues can be used as 
a diagnostic of the source of ionization and to 
distinguish the Hll-rcgion like objects from the 
Seyferts and LINERs. Most of the objects are 
found to be consistent with HII regions, e.g., 
from Liang et al. (2004) and exhibit low ion- 
ization level ([0III]/II;3 <3). From emission 
lines studies, the AGN fraction is quite low, 
~20 %. This is consistent with X-ray observa- 
tions of ISOCAM sources showing that AGNs 
contribute at most 17±6% of the total mid- 
infrared flux (Fadda et al. 2002). Assuming 
template SEDs typical of star-forming and star- 
burst galaxies, 15 fim fluxes can be converted 



into total infrared luminosities, Lir (between 8 
and 1000 ^m). About 75% of the galaxies dom- 
inated by the star formation are either LIRGs 
or ULIRGs. The remaining 25% are nearly 
equally distributed among either "starbursts" 
(101° < LiR < 10" Lq) or "normal" (Ljr < 
lO^" Lq) galaxies. The median luminosity is 
about 3x10" Lq. ULIRGs and LIRGs con- 
tribute to about 17% and 44% to the CIB at 15 
fj,m, respectively (Elbaz et al. 2002). This sug- 
gests that the star formation density at 2 <1 
is dominated by the abundant population of 
LIRGs. As will be shown later, this has im- 
portant consequences for the evolution of galax- 
ies. Because of large extinction in LIRGs and 
ULIRGs, the infrared data provide more robust 
SFR estimate than UV tracers. The extinction 
factor in LIRGs averages to Ay '^2.8 at z ~0.7 
(Flores et al. 2004). It is much higher than 
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that of the local star-forming galaxies for which 
the median is 0.86 (Kennicutt 1992). Assuming 
continuous burst of age 10-100 Myr, solar abun- 
dance, and a Salpeter initial mass function, the 
SFR can be derived from the infrared luminosi- 
ties (Kennicutt 1998): 

SFR{MQ/yr) = 1.71xlO"i"L/fl[8-lOOOMm](L0) 

(2) 

Thus typical LIRGs form stars at >20 
year^^. The median SFR for the 15 /im 
galaxies is about 50 M0 year"^, a substantial 
factor larger than that found for faint-optically 
selected galaxies in the same redshift range. 

The other fundamental parameter character- 
izing the sources of the peak of the infrared 
background is their stellar mass content that 
traces the integral of the past star forma- 
tion activity in the galaxies and is a natural 
complement to the instantaneous rate of star 
formation. The stellar masses can be obtained 
using spectral synthesis code modeling of 
the UV-optical-ncar infrared data or, more 
simply using the mass-to-luminosity ratio in 
the K-band. The derived stellar masses for the 
bulk of ISOCAM galaxies range from about 
10^° to 3x10"Mq, compared to 1.8x10" 
Mq for the Milky Way. As expected from 
the selection based on the LW3 flux limit - 
and thus on the SFR - masses do not show 
significant correlation with redshift (Frances- 
chini et al. 2003). An estimate of the time 
spent in the starburst state can be obtained by 
comparing the rate of ongoing star formation 
(SFR) with the total mass of already formed 
stars: t[years\ = M/SFR. Assuming a constant 
SFR, t is the timescale to double the stellar 
mass. For LIRGs at z >0.4, t ranges from 
0.1 to 1.1 Gy with a median of about 0.8 
Gyr (Franceschini et al. 2003; Hammer et al. 
2005). From z = 1 to z = 0.4 (i.e., 3.3 Gyr), 
this newly formed stellar mass in LIRGs corre- 
sponds to about 60% of the z — 1 total mass 
of intermediate mass galaxies. The LIRGs are 
shown to actively build up their metal content. 
In a detailed study, Liang et al. (2004) show 
that, on average, the metal abundance of 
LIRGs is less than half of the z ^ disks with 
comparable brightness. Expressed differently, 
at a given metal abundance, all distant LIRGs 
show much larger B luminosities than local 



disks. Assuming that LIRGs eventually evolve 
into the local massive disk galaxies, Liang et 
al. (2004) suggest that LIRGs form nearly half 
of their metals and stars since z ~1. 

Finally, morphological classification of dis- 
tant LIRGs is essential to understand their for- 
mation and evolution. Zheng et al. (2004) 
performed a detailed analysis of morphology, 
photometry, and color distribution of 36 (0.4< 
z <1.2) ISOCAM galaxies using HST images. 
Thirty-six percents of LIRGs are classified as 
disk galaxies with Hubble type from Sab to 
Sd; 25% show concentrated light distributions 
and are classified as Luminous Compact Galax- 
ies (LCGs); 22% display complex morphology 
and clumpy light distributions and are classi- 
fied as irregular galaxies; only 17% are major 
ongoing mergers showing multiple components 
and apparent tidal tails. This is clearly differ- 
ent from the local optical sample of Nakamura 
et al. (2004) in the same mass range in which 
27%, 70%, <2%, 3% and <2% are E/SO, spi- 
rals, LCGs, irregulars and major mergers re- 
spectively. Consequences for galaxy evolution 
will be given in Section ing For most compact 
LIRGs, the color maps reveal a central region 
strikingly bluer than the outer regions. These 
blue central regions have a similar size to that 
of bulges and a color comparable to that of star- 
forming regions. Because the bulge/central re- 
gion in local spiral is relatively red, such a blue 
core structure could imply that the galaxy was 
forming the bulge (consistent with Hammer et 
al. 2001). It should be noticed that they find 
all LIRGs distributed along a sequence that re- 
lates their central color to their compactness. 
This is expected if star formation occurs first in 
the center (bulge) and gradually migrate to the 
outskirts (disk), leading to redder colors of the 
central regions as the disk stars were forming. 

4.2. Cosmic Evolution 

Another remarkable property of the 15 fim 
sources is their extremely high rates of evolu- 
tion with redshift exceeding those measured for 
galaxies at other wavelengths and comparable 
to or larger than the evolution rates observed 
for quasars. Number counts at 15 /im show a 
prominent bump peaking at about 0.4 mJy. 
At the peak of the bump, the counts are one 
order of magnitude above the non-evolution 
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models. In fact, data require a combination 
of a (1+z)^ luminosity evolution and (1+z)^ 
density evolution for the starburst component 
at redshift lower than 0.9 to fit the strong 
evolution. Although it has not been possible 
with ISOCAM to probe in detail the evolution 
of the infrared luminosity function, Spitzer 
data at 24 /xm give for the first time tight 
constraints up to redshift 1.2 (Le Floc'h et al. 
2005; Perez- Gonzalez ct al. 2005). A strong 
evolution is noticeable and requires a shift of 
the characteristic luminosity L* by a factor 
(l+z)4-0±o.5^ Lc Floc'h ct ai. (2005) find that 
the LIRGs and ULIRGs become the dominant 
population contributing to the comoving in- 
frared energy density beyond z ^0.5-0.6 and 
represent 70% of the star-forming activity 
at z ~1. The comoving luminosity density 
produced by luminous infrared galaxies was 
more than 10 times larger at z ^1 than in the 
local Universe. For comparison, the B-band 
luminosity density was only three times larger 
at z = 1 than today. Such a large number 
density of LIRGs in the distant Universe 
could be caused by episodic and violent star- 
formation events, superimposed on relatively 
small levels of star formation activity. This 
idea has emerged in 1977 (Toomre 1977) and 
is fully developed in Hammer et al. (2005). 
These events can be associated to major 
changes in the galaxy morphologies. The 
rapid decline of the luminosity density from 
^; = 1 is only partially due to the decreasing 
frequency of major merger events. Bell et 
al. (2005) showed that the SFR density at 
z ~0.7 is dominated by morphologically normal 
spiral, E/SO and irregular galaxies (>70%), 
while clearly interacting galaxies account for 
<30%. The dominent driver of the decline is 
a strong decrease in SFR in morphologically 
undisturbed galaxies. This could be due, for 
example, to gas consumption or to the decrease 
of weak interactions with small satellites that 
could trigger the star formation through bars 
and spiral arms. 

Locally 0.5% of galaxies with Ly >10^^ Lq 
have SEDs typical of LIRGs. This changes dra- 
matically at higher redshift: in deep surveys, 
ISO detect about 15% of the Mb <-20 galaxies 
(LIRGs, Hammer et al. 2005) and Spitzer de- 



tect about 30% of field galaxies (Starbursts and 
LIRGs, Bell et al. 2005). Thus the two popu- 
lations (optical and LIR galaxies) overlap more 

at high z. 

4.3. Towards a Scenario of Recent Bulge and 
Disk Formation in Intermediate-Mass 
Spirals 

Because a significant part of recent star for- 
mation takes place in LIRGs, any overall picture 
of galaxy evolution requires a detailed panchro- 
matic study. Optical/spectral properties of 
LIRGs are similar to those of other galaxies and 
only infrared measurements are able to describe 
how the star formation is distributed between 
the different galaxy types. Thus a complete 
study has to link the star formation revealed in 
the infrared to the morphological changes seen 
in the optical. This has been done by Hammer 
et al. (2005) using HST, ISO, VLA and VLT 
observations of the CFRS. A detailed compari- 
son of the morphologies of distant (0.4<z<1.2) 
galaxies with the local galaxies shows the com- 
plete vanishing of the LCGs in the local Uni- 
verse (by a factor ~10) and the decrease of 
mergers and irregulars (by a factor ^4). Al- 
most all distant galaxies have much bluer cen- 
tral colors than local bulges, probably as a re- 
sult of active star formation in the Ikpc central 
region of most distant spirals. This supports 
a relatively recent formation of bulges in many 
present-day spirals. This simultaneous changes 
in galaxy morphologies and central colors of dis- 
tant galaxies together with the observed lower 
mctallicities (Liang et al. 2004) and overall 
higher star-formation rates at high z are the in- 
gredients for an updated scenario of bulge and 
disk formation in spirals. Hammer et al. (2005) 
propose three different phases of galaxy evo- 
lution: the mergers/interacting, the compact 
galaxy and finally the growth of disk phase. 
During the last 8 Gyrs, most luminous galax- 
ies are expected to experience a major merger 
that suppresses the disk as matter is falling to 
the mass barycenter. This phase is associated 
with short (1 Gyr) and strong peaks of star 
formations. Most of galaxies in this phase are 
LIRGs. Then, the compact phase corresponds 
to a decrease over 0.6-2 Gyr of the enhanced 
star formation due to merging. A significant 
fraction of stars form in bulges and additional 
occurrence of gas infall may subsequently wrap 
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around the bulge to form a new disk-like com- 
ponent. Finally, the star formation spreads over 
all the forming disk as fed by large amounts of 
gas infall. In this scenario, about half of the 
bulge stellar content was made earlier in their 
progenitors, before the last major phase of ac- 
cretion. More than a third of the present-day 
stellar mass is formed at z <1. This scenario is 
in very good agreement with the hydrodynami- 
cal numerical simulations of Scannapieco & Tis- 
sera (2003) in which mergers, through secular 
evolution and fusions, transform galaxies along 
the Hubble sequence by driving a morpholog- 
ical loop that might also depend on the prop- 
erties of the central potential wells, which are 
also affected by mergers. This very attractive 
scenario links in a simple way the distant and 
local galaxies; it will be confronted to the new 
panchromatic studies of Spitzer galaxies. Note 
that another possibility of buildup of dense cen- 
tral component in disk galaxies is internal secu- 
lar evolution, as reviewed by Kormendy & Ken- 
nicutt (2004). 

5. GALAXIES AT REDSHIFTS Z > 1.5 

Analysis of the CIB in the light of the ISO ob- 
servations shows that, as we go to wavelengths 
much longer than the emission peak, the CIB 
should be dominated by galaxies at higher red- 
shifts as illustrated in Figure 01 The comoving 
infrared production rate needed to fill the CIB 
around 1mm at a redshift centered around 2.5 
to 3 remains comparable to the one from galax- 
ies detected in the ISO CAM surveys and filling 
65% of the peak of the CIB. In this section we 
discuss the rapidly growing observational evi- 
dence that this picture is basically correct. The 
main source of these observations has been the 
SCUBA submillimeter observations at 850 fim 
and 450 /xm (see Blain et al 2002 for a re- 
view) and observations from the MAMBO in- 
strument on the IRAM 30-m telescope at 1.2 
mm (Greve et al. 2004). The negative K- 
correction becomes very effective at these wave- 
lengths, leading to an almost constant observed 
flux for galaxies of the same total infrared lumi- 
nosity between redshifts 1 and 5. More recently, 
the Spitzer observatory has produced a wealth 
of early observations showing that this observa- 
tory will contribute much to our understanding 
of infrared galaxies at z>1.5. 



5.1. Number Counts, Contribution to the CIB 

Blank-field deep surveys combined with map- 
ping of areas lensed by clusters lead to num- 
ber counts at 850 /im down to 0.5 mJy (e.g.. 
Small et al. 2002; Wang et al. 2004). At 
1.2 mm counts have been obtained down to 
2.5 mJy (e.g., Greve et al. 2004). The num- 
ber counts shapes at 850 /im and 1.2 mm are 
compatible with the assumption that they are 
made of the same population with a flux ratio 
-F85o/^'i20o = 2.5. For a typical ULIRG SED, 
a 5mJy source at 850 /xm has a luminosity of 
1O^^L0 at a redshift of about 2.5. The large 
fraction of the background resolved at 850 fj,m 
(see Section rO| has interesting consequences. 
It shows very directly that if the sources are at 
redshift larger than 1 (as confirmed by the red- 
shift surveys discussed below), the infrared lu- 
minosity of the sources that dominate the back- 
ground is larger than 10^^ Lq. This is a popu- 
lation with a very different infrared luminosity 
function than the local or even the 2: = 1 lu- 
minosity function. The link between this pop- 
ulation at high z, and what has been seen at 
2; '--^ 1 (as discussed in Section 0J will be done 
by Spitzer/MIPS observations at 24 /im. Fig- 
ure 21 shows that the building on the bulk of 
the CIB near its peak (at 150 /im) with red- 
shift is expected to be similar to the building 
of the 24 /im background when the history of 
the 15 and 70 fj,m CIB have larger contribu- 
tions from redshift-1 sources. The K-correction 
plots (Figure O show for 15 /im a hump at 
z = 1 associated with the coincidence of the 
6-9 fim aromatic features in the ISOCAM filter 
and a hump at the same redshift for the 24 /im 
MIPS filter associated with the 11-14 /im set of 
aromatic features in the MIPS filter. For the 
MIPS filter a second hump is visible at z ~2 
that corresponds to 6-9 /im features centered 
on the 24 /xm MIPS filter. ISOCAM galax- 
ies contribute to about 2/3 of the energy peak 
of the CIB. Following the previous considera- 
tions, it is easy to understand why the remain- 
ing fraction is likely to be made of sources in the 
redshift range 1.5-2.5. The presently detected 
SMGs with luminosity lO^^L© have an almost 
constant flux between redshift 1.7 and redshift 
2.5 at 24 /im (similar to the constant flux at 
850 /im between redshift 1 and 5). The MIPS 
24 /tm deep surveys (e.g., Papovich et al. 2004) 
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Fig. 7. — The observed histogram of the redshift distribution for the 55 radio-identified SCUBA galaxies {red 
histogram). Curves derived for a model of the radio/submillimeter galaxy populations (Chapman et al. 2003b; Lewis 
et al. 2005) are overplotted suggesting that the redshifts of the sources missed in the radio identification process 
lie mostly at redshifts ~3-5 between the radio and submillimeter model tracks. A sample of radio-selected QSOs is 
also overplotted {gray dotted line), revealing a remarkable similarity with the observed distribution for submillimeter 
galaxies. From Chapman et al. 2003a. 



reach a sensitivity of 50 /xJy and thus can de- 
tect all these galaxies when they are starburst- 
dominated. Considering the speed of the MIPS 
it is likely that 24 /im surveys will become the 
most efficient way to search for luminous star- 
burst galaxies up to z = 2.5 and up to 3 for the 
most luminous ones. 

5.2. Redshift Distribution and SEDs of the 
SMGs 

The first obvious question when investigat- 
ing the nature of the submillimeter galaxies 
(SMGs) is their redshift distribution. The 
rather low angular resolution of the submillime- 



ter and millimeter observations made identifi- 
cations with distant optical galaxies an almost 
impossible task without an intermediate identi- 
fication. This is provided by radio sources ob- 
served with the VLA with 10 times better an- 
gular resolution. The tight correlation between 
far-infrared luminosity and radio flux (Helou et 
al. 1985; Condon 1992) provides the needed 
link. This then allows us to get optical identifi- 
cations and redshift measurements using 10-m 
class telescopes. Confirmation of these identifi- 
cations can then be obtained through CO line 
observations with the millimeter interferome- 
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ters such as the Plateau de Bure interferome- 
ter. The redshift deduced from the optical lines 
is confirmed by the CO observations. So far, 
only a handful of cases have gone through this 
whole chain of observations (e.g., Genzel et al. 
2003; Greve et al. 2005; Neri et al. 2003), but 
a high success rate gives confidence in the first 
step of the identification process. The chain will 
also have to be applied to the tentative coun- 
terparts of radio-undetected SMGs that have 
been found using a certain combination of op- 
tical properties (Pope et al. 2005). 

The difficulty of making large, blind surveys 
at 850 fim at the required sensitivity has lead 
to an attempt to find distant SMGs through 
blind surveys at different wavelengths. Barger 
et al. (2000) have observed optically faint 
radio sources at submillimeter wavelength 
and demonstrated them to be, so far, the 
most efficient way to preselect targets for 
submillimeter observations and to get larger 
samples of potentially high redshift SMGs. As 
an example. Chapman et al. (2002) recovered 
at 850 fj,m 70% of the blind submillimeter 
survey sources. This contrasts the recovery 
rate of MAMBO sources, which is relatively 
low, --25% (Dannerbauer et al. 2004). It 
should be noticed that the radio preselection 
biases the sample against very high redshifts 
{z <3) because the radio flux at 1.4 GHz is 
below the detection limit of the VLA sur- 
veys used for this preselection. A model by 
Chapman et al. (2003b) and by Lewis et al. 
(2005) illustrate this effect very well (Figure 
I?!). A fraction of the submillimeter-selected 
sources are missed in such a process at z > 
3 (detectability in radio) and around z ~1.5 
(optical redshifts desert). The number of non 
identified submillimeter sources (around 30% 
for 5*850 > 3 mjy) is consistent with this 
model. Nevertheless the submillimeter-selected 
sources do not appear qualitatively different 
from the optically-faint-radio selected ones. 
Another bias is the effect of the dust effective 
temperature of the SMGs (Lewis et al. 2005). 
At a given total far-infrared luminosity, hotter 
sources have lower submillimeter fluxes if the 
radio/far- infrared correlation continues to hold. 
They could be missed in the submillimeter 
surveys (see the discussion in Chapman et al. 
2005). 



Chapman et al. (2003a) got spectroscopic 
redshifts of 55 sources obtained in this way. 
The redshift distribution for these sources is 
shown in Figure H (note that when this review 
was being edited. Chapman et al. (2005) pub- 
lish spectroscopic redshifts for 73 SMGs). This 
distribution peaks at z— 2.4 with a substantial 
tail up to z=4. In fact the redshift distribution 
can be represented by a Gaussian distribution 
centered on 2.4 and with a sigma of 0.65. 
Almost all SMGs are found in the redshift 
range 1.5 < z < 3. This redshift distribution is 
compared with that of the redshift distribution 
for a pure radio sample in Figure 13 The 
SMGs selected in the way described above is 
also shown to be very similar to the redshift 
distribution of the radio-selected QSOs. This 
observation is interesting in the context of high 
rate of AGN activity detected in SMGs. 
The determination of the SED of millime- 
ter/submillimeter galaxies remains an open 
question despite a lot of work in the last few 
years. The SCUBA and MAMBO data provide 
constraints on the flux and spectrum at long 
wavelengths; Spitzer observations constrain 
the near and mid-infrared. The far-infrared 
part of the SED remains the least precisely 
known. Low angular resolution makes 70 and 
160 fiia deep surveys confusion- limited at 3 
and 40 mJy (Dole et al. 2004b). These hmits 
are too high to complete the SED of the SMGs 
(see Figure EJ. Stacking sources will help 
to go deeper than the confusion limit when 
large samples of SMGs are available in MIPS 
cosmological surveys. A first attempt on a 
radio-selected sample lowered the limit down 
to 1.2 mJy at 70 //m (Frayer et al. 2004). 
They find a typical flux ratio I(70)/I(24)<7 
that they interpret to be low when compared 
with low-redshift starburst. However, such low 
ratios are typical of dusty starbursts placed at 
redshift greater than 1.5. It is thus likely that 
the lower colors are due to a redshift effect. 
Appleton et al. (2004) looked at the mid- and 
far-infrared fluxes from a purely radio-selected 
1.4 GHz fiJy sample of about 500 and 230 
sources at 24 and 70 fim, respectively. They 
show that the far-infrared to radio correlation 
that is constant out to z = 1 seems to be 
constant using 24 /im out to z = 2 but with a 
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larger dispersion due systematic variations in 
SED shape throughout the population. This 
provides positive evidence of the universality of 
the infrared/radio correlations out to redshifts 
of about 2. 

Blain ct al (2004a) have analyzed SEDs of in- 
frared galaxies assuming that the low-rcdshift 
radio/far- infrared correlation applies to SMGs. 
Under this reasonable assumption and using a 
model of long- wavelengths SEDs based on a sin- 
gle modific^d black body, they can choose a sin- 
gle parameter to built an SED that fits the long- 
wavelengths data and the radio/infrared lumi- 
nosity ratio. In their paper, this single param- 
eter is the temperature, but it could equally 
well be the long- wavelength emissivity, because 
they showed that this is degenerate with tem- 
perature. A split between two redshift popu- 
lations appear in their analysis. The high-^; 
galaxies selected by the submillimcter obser- 
vations are significantly colder that the low-z 
galaxies (Dune & Eales, 2001; Stanford et al. 
2000), IRAS or IRAS-radio selected. The dis- 
crepancy in part probably reflects selection ef- 
fects in the way these samples were obtained 
and may reflect the fact that SMGs and local 
infrared galaxies are distinct populations. It re- 
mains an open question what effect this has on 
the SED model. The main worry is that a single 
modified black body often does not fit ULIRGs 
SED when they are known at many frequencies. 
The SED is broader; the unavoidable temper- 
ature distribution of dust in infrared galaxies 
would affect such an analysis. In fact, the Stan- 
ford et al. (2000) sample does not agree well 
with the single-temperature SED, and this led 
Lagache et al. (2004) to take broader SEDs for 
their starburst galaxy templates. 

5.3. Nature of the SMGs 

Many LIRGs and ULIRGs at low redshifts 
have been identified with interacting or galaxy 
mergers. A substantial fraction show signs 
of AGN activity but it has been shown for 
the low-redshift LIRGs and ULIRGs that 
the starburst component dominates the en- 
ergy output (Genzel et al. 1998; Lutz et 
al. 1998). The sources used for the redshift 
distribution by Chapman et al. (2003a) have 
been imaged with the HST. Most of them 
are multi-component-distorted galaxy systems 



(Conselice et al. 2003; Small et al. 2004). They 
display irregular and frequently highly complex 
morphologies compared to optically selected 
galaxies at similar redshifts. They are often red 
galaxies with bluer companions, as expected 
for interacting, star-forming galaxies. They 
have higher concentrations, and more prevalent 
major-merger configurations than optically- 
selected galaxies at z ~2-3. Most strikingly, 
most of the SMGs are extraordinarily large 
and elongated relative to the field population 
regardless of optical magnitude (Chapman 
et al. 2003c). SMGs have large bolometric 
luminosities, ^ 10^^ — 10^'^ Lq, characteristic 
of ULIRGs. If the far-infrared emission arises 
from the star formation, the large luminosities 
translate to very high SFR >1000 M© year~^ 
Such high rates are sufficient to form the stellar 
population of a massive elliptical galaxy in 
only a few dynamical times, given a sufficient 
gas reservoir. SMGs are very massive systems 
with typical mass of l-2xlO^^L0 (Swinbank 
et al. 2004), comparable to the dynamical 
mass estimates from CO observations. Genzel 
et al. (2004; and more recently Greve et al. 
2005) have undertaken an ambitious program 
to study the nature of the SMGs in more 
details. They got CO spectra with the Plateau 
de Bure interferometer for 7 sources out of 
their sample of 12 for the CO 3-2 and 4-3 
transitions redshifted in the 3 mm atmospheric 
window. They provide optical identifications 
and redshifts. The detection of these sources 
at the proper redshift confirms the usefulness 
of identification with the help of the radio 
sources. The median redshift of this sample is 
2.4. In addition, one source was studied with 
the SPIFI instrument on the ESO/VLT. These 
observations are giving very interesting clues on 
the nature of the submillimeter galaxies. The 
gas masses obtained for these systems using 
CO luminosity /mass of gas determined from 
local ULIRGs is very large with a median of 
2.2x10^°Mq (10 times larger than in the Milky 
Way). Using the velocity dispersion, they could 
infer that the dynamical median mass of these 
systems is 13 times larger than in Lyman-break 
galaxies (LBGs) at the same redshift or 5 
times the mass of optically selected galaxies 
at this redshift. These SMGs with a flux at 
850 fxm larger than 5 mJy are not very rare 
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Fig. 8. — Comoving number densities of galaxies with baryonic masses > 10^^ A/0 as a function of redshift. The 
triangle and open squares show densities of massive stellar systems at 2 = and 2 ~ 1; The circle shows the density 
for massive SMGs at 2 ~ 2.7, with a factor of 7 correction for burst lifetime. Blue and red curves show the predictions 
of semianalytic models by the "Munich" and "Durham" groups, respectively. Dashed curves show the corresponding 
number densities of halos with available baryonic masses > 10" Mq. The two models use the same halo simulations 
but assume different f2b. From Genzel et al. (2004). 



and unusual objects, because they contribute 
to about 20% of the CIB at this frequency. 
Through multi wavelength observations, Genzel 
et al. (2004) get the stellar component in K 
band, and infer the star-formation rate and 
duration of the star-formation burst. They 
can then compare the number density of these 
massive systems with semiempirical models of 
galaxy formation. The very interesting result is 
that this number density is significantly larger 
than the predicted one, although the absolute 
numbers depends on a number of assumptions 
like the IMF. The comparison is shown in Fig- 
ure IHl Such massive systems at high redshift 
are not easy to understand in current cold 



dark matter hierarchical merger cosmogonies. 
However, one must keep in mind that bright 
SMGs (Sgso >5 mJy) that contribute 20% of 
the CIB may not be representative of the whole 
population. Gravitational lens magnification 
provides a rare opportunity to probe the nature 
of the distant sub-mJy SMGs. Kneib et al. 
(2005) study the property of one SMG with 
an 850 /im flux S85o=0.8 mJy at a redshift of 
z = 2.5. This galaxy is much less luminous and 
massive than other high-z SMGs. It resembles 
to similarly luminous dusty starbursts resulting 
from lower-mass mergers in the local Universe. 

In order to link the different population of 
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high-redshift objects, several LBGs at redshift 
between 2.5 and 4.5 have been targeted at 
850 /im. The Lyman-break technique (Steidel 
et al. 1996) detects the rest- frame 91.2 nm neu- 
tral hydrogen absorption break in the SED of a 
galaxy as it passes through several broad-band 
filters. LBGs are the largest sample of spec- 
troscopically confirmed high-redshift galaxies. 
Observing LBGs in the submillimeter is an im- 
portant goal, because it would investigate the 
link, if any, between the two populations. How- 
ever, the rather low success rate of submillime- 
ter counterpart of LBGs (e.g.. Chapman et al. 
2000; Webb et al. 2003) argues against a large 
overlap of the two populations. 

5.4. Spitzer 24 fJ-m Sources 

A potential new way to find high-z LIRGs 
and ULIRGs appeared recently with the launch 
of the Spitzer observatory. Particularly suited 
to this goal is the 24 /im channel of the MIPS 
instrument. The confusion levels in the 70 and 
160 /im prevent detection a significant number 
of high-redshift objects, and the IRAC 3.6 to 
8 /im at high redshift probes mostly the old 
stellar component that is much weaker than the 
dust emission in starburst galaxies. At the time 
of writing, the observations are under way, and 
only a few results are available. Le Floc'h et 
al. (2004) give the first hint on the 24 /im se- 
lected galaxies. They couple deep 24 /im obser- 
vations in the Lockman hole and extended groth 
strip with optical and near-infrared data to get 
both identification and redshift (either spectro- 
scopic or photometric). They find a clear class 
of galaxies with redshift 1< z <2.5 and with 
luminosities greater than ^5x10^^ L0 (see also 
Lonsdale et al. 2004). These galaxies are rather 
red and massive with M>2xl0^° Mq (Caputi 
et al. 2005). Massive star- forming galaxies re- 
vealed at 2 < z < 3 by the 24 /im deep sur- 
veys are characterized by very high star for- 
mation rates - SFR >500 Mq year'^ They 
are able to construct a mass of ~10^^ Mq 
in a burst lifetime (~0.1 Gyr). The 24 /im 
galaxy population also comprises sources with 
intermediate luminosities (10^" < Lj^ <10^^ 
Lq) and low to intermediate assembled stel- 
lar masses (10^ <M<10" Mq) a,t z < 0.8. 
At low redshifts, however, massive galaxies are 
also present, but appear to be building their 
stars quiescently in long timescales (Caputi et 



al. 2005). At these redshifts, the efficiency 
of the burst-like mode is limited to low mass 
M<10^° Mq galaxies. These results support a 
scenario where star-formation activity is differ- 
ential with assembled stellar mass and redshift, 
and proceed very efficiently in massive galaxies 
(Caputi et al. 2005). 

In the Lockman Hole, only one galaxy is 
associated with an X-ray source. This suggests 
that these galaxies are mostly dominating by 
star formation, consistent with the findings of 
Alonso-Herrero et al.(2004) and Caputi et al. 
(2005). This is also suggested by SEDs that 
are best fitted by PAH features rather than 
by strongly rising, AGN-type continua (Elbaz 
et al. 2005). The selected sources exhibit a 
rather wide range of MIPS to IRAC flux ratio 
and optical/near- infrared shapes, suggesting 
a possibly large diversity in the properties of 
infrared galaxies at high redshift as noticed 
by Yan et al. (2004b). Based on these first 
analyzes, together with the interpretation 
of the number counts (e.g., Lagache et al. 
2004), it is clear that the 24 /im observations 
will provide the sample to unambiguously 
characterize the infrared galaxies up to z ~ 2.5. 
They should fill the gap between the ISO- and 
SCUBA-selected galaxies. 

Several 24 /xm observations have been con- 
ducted on selected ERO and SCUBA and 
MAMBO samples. To our knowledge, LBGs 
have not been observed at long wavelengths. 
The MAMBO/SCUBA selected galaxies in the 
Lockman hole with radio identification have 
been observed by Spitzer and most of them 
detected between 3.6 and 24 /im. This allows 
to get an average SED for these (Egami et 
al. 2004; Ivison et al. 2004; see Figure ^ 
Spitzer deep surveys at 24 /im and shallow 
surveys like the SWIRE legacy (Lonsdale et 
al. 2004) can easily detect them and are thus 
a promising new way to find this class of 
high-z infrared galaxy. Nevertheless, the Early 
Release Observations from Spitzer have been 
used to extract their submillimeter flux from 
a stacking analysis of SCUBA observations in 
the Lockman hole (Serjeant et al. 2004). In 
this field, seven SMGs were already known 
and others were identified by further analysis. 
For the bulk of the 24 /im sources a marginal 




Fig. 9.— Rest-frame SED of 15 SMGs (assuming a redshift of 3) with MAMBO and/or SCUBA, Spitzer/IRAC 
and Spitzer/MTPS 24 fim measurements. Purple diamonds are the galaxies 208, 119, 115, 48, 44 (Frayer et al. 2004), 
LE850.4, LE850.35 (Egami et al. 2004), and MMJ105201, MMJ105155, MMJ105203, MMJ105216, MMJ105148, 
MMJ105157, MMJ105207 MMJ105203 (Ivison et al. 2004). Overplotted are the SEDs of M82, normaUzed at 850 
(from Chanial 2003), and the SED template of the Lagache et al. (2004) model, for L=10^'' Lq and a redshift of 3 
(no normalization has been applied). Note that this sample of SMGs has a ratio dust /stellar component higher than 
the template or M82. 



detection is found with an 585o/'S'24 ratio 
(1/20) much lower than that observed for 
SMGs. This clearly shows that the SMGs 
are only a fraction of the 24 fim sources, as 
expected. An interesting challenge is to find if 
Spitzer color criteria can be found to extract 
preferentially SMGs, i.e., the galaxies that 
account for most of the GIB near 1 mm. The 
SED in the thermal infrared appears quite 
variable for LIRGs and ULIRGs making this 
difficult (e.g., Armus et al. 2004). 

Extremely Red Objects (EROs) are usu- 
ally selected based on their red colors: {R — 
Ks) >5.3 mag or (/ — Kg) >4 mag. This color 
selection should include early-type galaxies at 
z However, the color selections are also sen- 



sitive to dust-reddened, star-forming systems. 
Up to now, it remains unclear what fraction of 
EROs are truly dust-obscured galaxies. Differ- 
ent scenarios of galaxy formation predict very 
different formation epochs for such galaxies. It 
is thus interesting to characterize these galax- 
ies, in particular whether they belong to the 
early-type or dusty star-forming class of ob- 
jects. Spitzer/MlVS 24 /im observations offer 
the first opportunity to address this issue be- 
cause 24 /im observations can clearly discrim- 
inate between the two populations. In the Nl 
field, Yan et al. (2004a) suggest that about 50% 
of EROs are infrared luminous, dusty starbursts 
at z > 1 (in a similar study, Wilson et al. (2004) 
show that at least 11% of 0.6 < z < 1.3 EROs 
and at least 22% of z > 1.3 EROs are dusty 
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star- forming galaxies). Their mean 24 /xm flux 
corresponds to infrared luminosities of about 
SxlO^'^ and 10^^ L0 at z^l and z ^1.5, respec- 
tively. They are massive galaxies with lower 
limit M>5xl09 to 2 10'^° Mq. The fraction 
of EROs likely to be AGN is small; about 15%. 
The link between the two classes of EROs could 
be that starburst EROs are experiencing, at 
z > 1, violent transformations to become mas- 
sive carly-typc galaxies. 

5.5. ULIRGs and Active Galactic Nuclei at 
High Redshifts 

SMGs are massive ULIRGs at high redshift. 
One of the key question discussed above for the 
~1 galaxies is to distinguish whether star- 
burst or AGN activity powers the dust heating 
and associated infrared emission. The presence 
of an AGN in galaxies can be investigated us- 
ing optical/near- infrared, emission line diagnos- 
tics and/or X-ray observations. But the iden- 
tification of the presence of an AGN does not 
mean that it is the dominant source of the far- 
infrared emission. Alexander et al. (2003; see 
also Almaini et al. 2003) use Chandra obser- 
vations of the CDF-N to constrain the X-ray 
properties of 10 bright SMGs. Half of the sam- 
ple has flat X-ray spectral slopes and luminous 
X-ray emission, suggesting obscured AGN ac- 
tivity. However, a comparison of the AGN- 
classified sources to the well-studied, heavily 
obscured AGN NGC 6240 suggests that the 
AGN contributes on average a negligible frac- 
tion (about 1.4%) of the submillimeter emis- 
sion. For the MAMBO sources, similar results 
are found: only one out of the nine MAMBO 
sources studied by Ivison et al. (2004) has an 
X-ray counterpart. It has, as expected from low 
redshift ULIRGs observations (e.g., Rigopoulou 
et al. 1999), a different mid- infrared SED than 
the starburst dominated sources. About 75% 
of their sample has rest-frame mid-infrared to 
far-infrared SED commensurate with obscured 
starburst. Swinbank et al. (2004), using AGN 
indicators provided by near-infrared spectra, es- 
timate that AGNs are present in at least 40% 
of the galaxies in their sample of 30 SMGs. 
Emission-line diagnostics suggest that star for- 
mation is the dominant power source. However, 
the composite spectrum for the galaxies that 
individually show no signs of an AGN in their 
near-infrared spectra appears to show an under- 



lying broad Hq line. This suggests that even 
these galaxies may host a low-luminosity AGN 
that is undetectable in the individual spectra. 
All these studies tend to show that starburst 
activity is the dominant source of power of dust 
emission in the far-infrared. Still, it is rather 
difficult to estimate the true "contamination" 
by the AGN. To go deeper, Chapman et al. 
(2004) tried an original approach. They ob- 
serve a sample of identified SMGs at high an- 
gular resolution in the radio and use the radio 
emission as a proxy for the far-infrared emis- 
sion. This assumption is based on the well- 
known very tight far-infrared/radio correlation 
mentioned above. If detected, an extended ra- 
dio (and thus far-infrared) component is likely 
to arise from the star formation. The detection 
of extended emission requires sub-arcsec reso- 
lution to map emission on kpc-scales. These 
are accessible by radio interferometry (they are 
well beyond far-infrared and submillimeter fa- 
cilities capabilities). They find that for 70% 
of the SMC sample, the MERLIN/VLA radio 
exhibits resolved radio emission which mirrors 
the general form of the rest frame UV morphol- 
ogy seen by HST. The galaxies are extended on 
scales of about 10 kpc. They interpret this as 
a strong support for the hypothesis that radio 
emission traces spatially extended massive star 
formation within these galaxies. This is clearly 
different from what is seen in local ULIRGs 
where the far-infrared/radio emission is con- 
centrated in the compact nuclear region with 
an extend less than 1 kpc. In the remaining 
30% of the SMG sample, the radio emission is 
more compact (essentially unresolved). This is 
a signature of either a compact nuclear star- 
burst and/or an AGN. 

In conclusion, the exact fraction of distant 
submillimeter and millimeter galaxies contain- 
ing an energetically dominant AGN is difficult 
to extract from observations. However, even in 
the systems containing an unambiguously pow- 
erfully AGN, the far- infrared emission seems 
to be powered by the star formation. Surpris- 
ingly, this seems to be also the case in distant 
QSOs. Recently, Beelen (2004) has shown that 
the far-infrared and blue luminosities from the 
host galaxies of distant radio-quiet QSOs, are 
slightly correlated. The far-infrared and ra- 
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Fig. 10. — Comoving correlation length of the SMGs 
{triangle) in contrast to other populations of low- and 
high-redshift galaxies (see the summary in Overzier et 
al. 2003). The horizontal error bar on the SMG point 
spans the range of redshifts over which SMG associa- 
tions are found. The solid line shows a representative 
model for the evolution of a certain overdensity. The 
dashed lines show the expected correlation length of 
dark matter halos as a function of mass and redshift. 
From Blain et al. (2004a). 



dio emission of these quasars follow the radio- 
infrared correlation observed in local ULIRGs 
(Yun et al. 2001), providing a first indica- 
tion that the dust is predominantly heated by 
the star-formation activity rather than by the 
AGN. Moreover, the non-linearity between the 
far-infrared and blue luminosities is also an in- 
dication that the heating mechanism of the dust 
is not directly linked to the AGN. However, 
the presence of this correlation could suggest 
a causal connection between the formation of 
stars in the host galaxy and the activity of the 
central super massive black hole. This connec- 
tion has been successfully modeled by Granato 
et al. (2004). 

Finally Houck et al. (2005) and Yan et al. 
(2005) demonstrate the potential of using mid- 
infrared spectroscopy, especially the aromatic 
and silicate features produced by dust grains 
to directly probe distant L~10^'^ Lq ULIRGs 
at z ~ 2. Spitzer /IRS observations provide 



a unique and direct access to high-z ULIRG 
physical properties. It will definitively open the 
route toward a complete census of the distant 
infrared-luminous Universe. A first study on 
two distant SMGs using Spitzer/IRS by Lutz et 
al. (2005) finds for one SMG an equal contribu- 
tion from star formation and AGN. The second 
galaxy is dominated by star formation. 

6. CLUSTERING 

Measuring clustering gives information about 
the distribution of galaxies with respect to 
the dark matter. The strength of clustering is 
correlated with the mass of extended halos that 
host luminous galaxies. At redshift lower than 
1.3, large-scale structures have been mapped by 
the DEEP2, 2dF and SDSS surveys (Coil et al. 
2004; Peacock et al. 2001; Doroshkevich et al. 
2004). At higher redshifts, correlation lengths 
ro have been measured for different galaxy 
populations such as the LBGs and EROs. 
Values of about ro=4 h~^ Mpc (Porciani & 
Giavahsco 2002) and ro = ll h^^ Mpc (Daddi et 
al. 2000; Firth et al. 2002; Roche et al. 2003) 
are reported with mean redshift of about 3 and 
1-1.5 for the LBGs and EROs, respectively. 
Getting information on the clustering of the 
infrared galaxies is essential to see how they 
relate to the other galaxy populations and to 
understand their formation process. As an 
example, one of the key questions is to see if 
the most massive SMGs are associated with 
the most massive dark matter halo. 

Up to now, very little information is available 
on the clustering of infrared galaxies. A clus- 
tering signal can tentatively be measured, al- 
though the small number of objects in the deep 
narrow-pencil ISOCAM and SCUBA beams 
prevents accurate measurements of the auto- 
correlation function. At z ~1, Moy & Elbaz 
(2005) study the large-scale and close environ- 
ment of ISOCAM galaxies. They find that 
infrared galaxies are more strongly clustered 
than optical galaxies. Eighty percents of ISO- 
CAM galaxies are found preferentially in red- 
shift peaks, versus 68% for the optically selected 
galaxies. Moy & Elbaz (2005) find indirect ev- 
idences that the triggering mechanism of dusty 
starbursts is small-scale 100 kpc) galaxy- 
galaxy interactions. Such interactions do not 
lead to major mergers most of the time but are 
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more likely simple fly-by, tidal interactions or 
minor mergers. At much higher redshift, there 
are some indirect evidences of strong cluster- 
ing of SMGs when compared to other classes 
of high-redshift galaxies. For example, Chap- 
man et al. (2001) identified SMGs within the 
most overdense structure of LBGs at z ~3.1 
(Steidel et al. 2000). De Breuck et al. (2004) 
detected an overdensity of MAMBO sources 
likely at z ~4.1 in a proto-cluster containing 
also overdensities of Lya emitters and LBGs. 
The first three-dimensional quantitative mea- 
surement of the clustering strength of SMGs has 
been made by Blain et al. (2004a). They find 
spectroscopic evidence for clustering. Using 73 
spectroscopically identified galaxies, they find 
a surprisingly large number of "associations" 
with redshifts separated by less than 1200 km 
s~^. They provide tentative evidence for strong 
clustering of SMGs at z ~2-3 with a corre- 
lation length of '^(6.9±2.1) h~^ Mpc using a 
simple pair-counting approach appropriate for 
the small and sparse SMC sample. This cor- 
relation length appears to be somewhat larger 
than that for both LBG and QSO galaxies at 
comparable redshifts. It is thus unlikely that 
the SMGs form a simple evolutionary sequence 
with either population. On the contrary, the 
correlation length could be consistent with a 
form of evolution that subsequently matches 
the large comoving correlation length typical of 
evolved EROs at z ~1 and of clusters of galax- 
ies at z = (Figure EB- From this figure, 
we see that the correlation function of SMGs 
appears to be consistent with the hypothesis 
that they are associated with the most mas- 
sive dark matter halos at high redshifts. These 
are more massive than the host halos of LBGs 
and QSOs at comparable redshift. These pre- 
liminary conclusions have to be investigated in 
more detail. In particular, the SMC masses 
inferred from clustering measurements have to 
be compared to the dynamical masses derived 
from millimeter wave CO spectroscopy. To go 
deeper in understanding the clustering prop- 
erties of infrared galaxies, Spitzer 24 /im sur- 
veys will be best suited. With a large num- 
ber of sources detected in quite large surveys, 
Spitzer will unambiguously constrain the clus- 
tering of infrared galaxies from redshift 0.5 to 
2.5. Conjointly, the physics of galaxy clustering 



can be probed by the GIB fluctuation analysis 
(Knox et al. 2001). CIB fluctuations measure, 
on large angular scales, the linear clustering 
bias in dark-matter halo and, at small angular 
scales, the nonlinear clustering within a dark- 
matter halo (Cooray & Sheth 2002). They thus 
probe both the dark-matter halo mass scale and 
the physics governing the formation of infrared 
galaxies within a halo. Promising attempts are 
underway using the 170 ^m ISO FIRBACK 
fields (Lagache et al., in preparation). 

7. FINDING ULIRGS AT Z >3: CONFUSION 

In the near future, when a proper census 
of ULIRGs up to z ~3 will have been carried 
out, the fraction of the CIB at ~1 mm not 
accounted for should give an indication of the 
contribution from sources at larger redshifts. 
Deep surveys at '^1-2 mm are the only obvious 
tool to find most of these sources. However, 
the limiting factor of the surveys is not only 
detector sensitivity or photon noise but also 
confusion. We concentrate in this section only 
on extragalactic sources confusion. Source 
confusion in the far-infrared is illustrated in 
Figure ^2 

Predicting or measuring confusion depends 
on the scientific goal of the measurement (Helou 
& Beichman 1990; Dole et al. 2003; La- 
gache et al. 2003). Performing an unbiased 
far-infrared or submillimeter survey and get- 
ting a complete sample has different require- 
ments than following-up in the far-infrared an 
already known near-infrared source to get an 
SED and/or a photometric redshift. In the for- 
mer case, one has to tightly control the statis- 
tical properties of the whole sample; in the lat- 
ter case, completeness is irrelevant, and even 
a low photometric accuracy is adequate. We 
thus favor the use of a term like "unbiased con- 
fusion" for the former case. New techniques 
are being developed to use a priori informa- 
tion at shorter wavelength (e.g., 8 /im with 
Spitzer/IRAC and 24 ^m with MIPS) to infer 
some statistical properties (such as source den- 
sity or SED) of sources at longer wavelength 
(e.g., 24 or 160 /im, respectively), and thus to 
beat unbiased confusion. Predicting unbiased 
confusion (for instance Condon 1974; Frances- 
chini et al. 1989; Helou & Beichman 1990; Dole 
et al. 2003; Lagache et al. 2003; Takeuchi et 
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Fig. 11. — Effects of confusion in the far-infrared. Observation of a source in the ELAIS-Nl/FIRBACK field in a 
400x400 square arcsec box. All images have been resampled to 5 arcsee per pixel, which oversamples the far-infrared 
maps but undersamples the mid-infrared map. From left to right: 373 mjy ISOPHOT 170 fim source with about 128s 
of integration (FIRBACK survey, Dole et al. 2001); labels indicate the 5 arcsec pixels; Spitzer/MTPS 160 ^im with 
about 16s of integration (SWIRE survey, Lonsdale et al. 2004); MIPS 70 fiui with about 80s of integration (SWIRE); 
MIPS 24 fim with about 160s of integration (SWIRE). Notice (1) the ISO 170 fim source is marginally resolved with 
MIPS 160, and is unambiguously resolved at 70 fira and 24 fira; (2) the two fainter MIPS 160 fim resolved sources 
(bottom left) create fluctuations in the ISO 170 fim map that produce the confusion noise when the resolution is 
limited. 



al. 2004; Negrello et al. 2004) requires the 
knowledge of at least the number-count dis- 
tribution of the galaxies. In practice, models 
(validated at some point by observations) are 
used. Because the shape of the counts in a 
log{N) — log{Si,) diagram varies with the flux 
density Sn, the fluctuation level of faint sources 
below Si, also vary. This fluctuation level gives 
an estimate of the unbiased confusion using a 
photometric criterion (Lagache & Puget 2000; 
Dole et al. 2003; Lagache et al. 2003). At very 
faint fluxes, when the background is almost re- 
solved, the photometric criterion will obviously 
give a very small value for the unbiased con- 
fusion level, but the observations will be lim- 
ited by the confusion due to the high density 
of faint resolved sources. Thus, another crite- 
rion, the source density criterion for unbiased 
confusion (SDC, Dole et al. 2003; Dole et al. 
2004b), needs to be computed and compared to 
the photometric criterion. 

In the infrared and submillimcter range be- 
low 300 /im, the unbiased confusion is in gen- 
eral better predicted by the source density cri- 
terion for current and future facilities, because 
the angular resolution has improved (e.g. from 
ISO to Spitzer). At longer wavelengths, the 
photometric criterion is more useful. We use 
the model of Lagache et al. (2004) to predict 
unbiased confusion. Fifteen-meter submillime- 



ter telescopes are limited by the confusion at 
1.2 mJy at 850 /xm and at 0.5 mJy at 1.2 mm. 
Ongoing surveys (Small et al. 2002; Greve et al. 
2004) already reach or are about to reach these 
levels. If we want to resolve about 80% of the 
GIB, which corresponds to 56 /xJy at 850 /zm 
and 20 /zJy at 1.3 mm, one would need a ~90- 
m telescope at 850 /im and a '^150-m telescope 
at 1.4 mm. Future facilities for infrared and 
submillimeter observations include far-infrared 
space observatories such as Herschel, SPIGA, 
and SAFIR, survey missions like ASTRO-F and 
Planck^ a larger near-infrared and mid-infrared 
observatory, JWST, and a ground-based sub- 
millimeter interferometer, ALMA. In order to 
detect LIRGs a.t z ~ 3, experiments operating 
at 15, 24, 70, 160, 450, 850, 1380 fim respec- 
tively should reach a sensitivity of 1, 1, 8, 50, 
200, 70, 30 fiJy, respectively (see Figure [T^ . 
The spectral window around 450 /xm seems the 
most effective to reach these galaxies. This con- 
straint is somewhat relaxed if one wants to de- 
tect ULIRGs at z ~ 3, in which case the re- 
quired sensitivities are multiplied by about 10. 
In the near future, ASTRO-F, Herschel/SFIRE 
and Planck will be mostly limited by confusion. 
At long wavelengths, to probe most of the GIB 
source population and to detect enough early 
mergers made by building blocks not yet af- 
fected by star formation and evolution, large 
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Fig. 12. — Sensitivity to the bolometric luminosity (and star-formation rate, assuming star-forming galaxies) of 
hypothetical surveys designed to detect LIRGs at z ~ 3. The required sensitivities are: at 15 lira: Su > l^Jy; at 
24 /tm: > l/^Jy; at 70 ^im: & > 8/xJy; at 160 fim: Si, > SO^tJy; at 450 ^ira: S„ > 200/^Jy; at 850 ^tm: > 70/xJy; 
at 1.3 mm: Si, > 30/iJy. This plot makes use of the Lagache et al. (2004) model (see appendix). 



extragalactic surveys will have to be conducted 
with ALMA. These surveys will take a substan- 
tial fraction of the time (Lagache et al. 2003). 
As an example, mapping one square degree at 
1.3 mm at the 5a sensitivity of 1 mJy - ~50% 
of the CIB is resolved - takes 138 days without 
including overheads. 

8. COMPARISON OF MODELS WITH 
OBSERVATIONS 

One of the striking result of the deep surveys 
concerns the evolution of the infrared and 
submillimeter galaxy population. The source 
counts are high when compared to no evolution, 
or moderate, evolution models^ for infrared 
galaxies. Classical semianalytical models of 
galaxy formation predicts neither the large 

^ 'No-evolution': the co-moving luminosity function 
remains equal to the local one at all redshifts 



numbers of infrared galaxies nor their very 
strong evolution, revealing a serious gap in 
our understanding of galaxy formation and 
evolution. Very recently, several empirical ap- 
proaches have been proposed to model the high 
evolution of the infrared output with redshift 
(e.g.. Chary & Elbaz 2001; Franceschini et al. 
2001; Rowan-Robinson 2001; Takeuchi et al. 
2001; Xu et al. 2001; Lagache et al. 2004) 
that fit source counts, redshift distributions 
and CIB intensities and fluctuations, although 
often not all of them. All these models, 
however, agree on a general trend - i.e., the 
luminosity function must change dramatically 
with redshift, with a rapid evolution of the 
high-luminosity sources (L>2 10^"'^ L0) from 
z = to z = 1, which then stay rather constant 
up to redshift 3 or more. The evolution of the 
infrared luminosity function may be linked to a 
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bimodal star-formation process, one associated 
with the quiescent and passive phase of the 
galaxy evolution and one associated with 
the starburst phase, triggered by merging 
and interactions. The latter dominates the 
infrared and submillimcter energy density of 
the Universe at high z. Consistently, cold dark 
matter N-body simulations show that halo 
merger rates increase with redshift as (1+z)™ 
with 2.5< m <3.5 (Gottlober et al. 2001). 
Observations, however, give m values between 
and 4 (Lc Fcvrc ct al. 2000; Consclicc ct 
al. 2003; Bundy ct al. 2004; Lin et al. 2004). 
The spread is due to dilferent selection effects, 
detection techniques, pair criteria and sample 
variance. It is therefore not easy to reconcile 
the different observational results. Moreover, 
comparisons with models are very difficult 
because definitions of merger rates may not be 
consistent. Merger rates can also depend on 
halo masses. As a consequence, the timescale 
of the merger phase is difficult to estimate. 
Peaks of star formation produced by mergers 
in hydrodynamical models (e.g., Scannapieco 
& Tissera 2003) has a duration of several hun- 
dreds of million years. This is consistent with 
what is observed. ULIRGs emit more than half 
of their bolometric luminosity from a starburst 
of age 10^-10* years (Genzel et al. 1998). 
LIRGs build up their stellar mass in a typical 
timescale of about 0.1 Gyr (Pranceschini et al. 
2003). These timcscalcs arc also supported by 
Marcillac et al. (2005) who performed Monte 
Carlo simulations using synthetic spectra based 
on the models of Bruzual and Chariot (2004) 
to derive the past star-formation history of 22 
LIRGs. They found that LIRGs experience a 
major event of star formation in their lifetime 
that produce about 10% of their stellar mass 
within 0.1 Gyr. How many such episodes of 
violent star formation does a typical galaxy 
experience? Assuming a timescale of 0.1 Gyr, 
Hammer et al. (2005) estimate the number of 
episodes per galaxies as about 5 from z = \ 
to 2: = 0.4. These episodic bursts naturally 
explain the high fraction of LIRGs in the 
distant Universe. 

Models that are more sophisticated than em- 
pirical approaches attempt to follow the physics 
of galaxy formation in greater detail (e.g., 



Guiderdoni et al. 1998; Hatton et al. 2003; 
Granato et al. 2004; Silva et al. 2005). In 
semianalytical models, the collapse of pertur- 
bations is described by the classical top-hat 
model under the assumptions of homogeneity 
and sphericity. The mass distribution of col- 
lapsed halos is computed from the so-called 
peaks formalism developed by Bardeen et al. 
(1986). Then dissipative collapse and cooling 
are introduced, with the usual "overcooling" 
problem that can partly be solved by introduc- 
ing stellar feedback. Star-formation processes 
are deduced from the gas content and the dy- 
namical timescale of the galaxies. Finally spec- 
trophotometric evolution is used to compute the 
age dependence of the gas content, the spec- 
tra of the stellar populations and the mass-to- 
luminosity ratios. To make specific predictions 
for the infrared galaxies, these models must in- 
clude an important additional feature: absorp- 
tion of the UV/optical radiation and emission 
by the dust grains. Very often two modes of star 
formations are considered; a quiescent mode 
and a burst mode in which the star formation 
timescales are much shorter. This burst mode 
is triggered by galaxy mergers and is absolutely 
required by the infrared to submillimcter ob- 
servations. There are some indications that to 
reproduce the submillimcter galaxy counts, a 
dramatic change of the IMF is required. A top- 
heavy IMF, in particular, increases the produc- 
tion of dust that is essential for boosting the lu- 
minosity of galaxies in the submillimcter. Using 
an IMF of the form dN/d\nm oc with x=0 
for the burst mode, Baugh et al. (2005) were 
able to reproduce not only the submillimcter 
observations but also the properties of Lyman- 
Break galaxies. They predict that the SMGs 
reside in the more massive halos in place at 
z = 2 and therefore that they are more strongly 
clustered than dark matter at this epoch. This 
is consistent with tentative observational con- 
straints (Blain et al. 2004a). There are several 
observational "indications" of massive stars (> 
100 Mq) in nearby starburst templates. Wolf- 
Rayet stars^ have been detected in a large num- 

2 Wolf-Rayet stars are hot (25,000 to 50,000 K), mas- 
sive (> 25M0), luminous stars with a high rate of mass 
loss. The Wolf-Rayet phase appears in an advanced 
stage of evolution. They are believed to be O stars that 
have lost their hydrogen envelopes, leading their helium 
cores exposed. Wolf-Rayet stars are often in a binary 
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ber of galaxies undergoing intense bursts of star 
formation (e.g., Gonzalez-Delgado et al. 1997; 
Pindao ct al. 2002). However, it remains diffi- 
cult to measure the IMF at high mass because 
of aging effects that can mimic real upper-mass 
IMF cutoff (the highest massive stars have very 
short lifetimes). 

In conclusion, the hierarchical galaxy formation 
paradigm is very successful in its description of 
large-scale structure formation and evolution. 
The next important step will be to test this 
picture to explain not only the number densities 
but also the mass assembly and particularly the 
mass of the SMGs. First mass measurements of 
SMGs galaxies seem to show that a very flat 
IMF cannot by itself explain the mass assem- 
bly of the baryonic matter at high z (Genzel 
et al. 2004). Hierarchical clustering underpre- 
dicts the high-z volume densities of these mas- 
sive galaxies. More work needs to be done to 
test the baryonic mass assembly in the hierar- 
chical paradigm. Both observational and model 
estimates are still very uncertain, with the for- 
mer depending on large lifetime corrections and 
small samples and the latter on ad hoc input 
recipes for feedback and star formation. 

9. CONCLUSION AND OPEN QUESTIONS 

A number of conclusions are now clear from 
the analysis of the identified sources in the GIB: 

• The comoving energy produced in the 
past that makes up the GIB at different 
wavelengths is more uniform that what is 
suggested by its spectral energy distribu- 
tion. This is due to the fact that the GIB 
at long wavelengths (A > 400/Ltm) is dom- 
inated by emission from the peak of the 
SED of galaxies at high z. More quanti- 
tatively, the ISOGAM surveys reveal that 
about two-thirds of the GIB emission at 
A ~150 /zm is generated by LIRGs at 
z ~ 0.7. At 850 /zm, more than half of the 
submillimeter GIB is generated by SMGs. 
The brightest SMGs (Sgso >3 mJy, -30% 
of the GIB) are ULIRGs at a median red- 
shift of 2.2. The energy density at 150 /im, 
which is ~20-25 times larger than the en- 
ergy density at 850 /im requires a comov- 
ing energy production rate zk. z = 0.7 

system, and are deemed, within a few million years, to 
explode as type lb or la supernovae. 



roughly 10 times the energy production 
rate at .2 = 2.2. 

• The evolution exhibited by LIRGs and 

ULIRGs is much faster than for optically 
selected galaxies. The ratio of infrared to 
optical, volume-averaged output of galax- 
ies increases rapidly with increasing red- 
shift. 

• Luminosity function evolution is such that 
the power output is dominated by LIRGs 
at 2; ~ 0.7 and ULIRGs at z ~ 2.5. 

• The energy output of GIB sources is dom- 
inated by starburst activity. 

• AGN activity is very common in the most 
luminous of these galaxies even though 
this activity does not dominate the en- 
ergy output. The rate and fraction of the 
energy produced increase with the lumi- 
nosity. 

• LIRGs at 2; ~ 0.7 are dominated by inter- 
acting massive late-type galaxies. Major 
mergers become dominant in ULIRGs at 
z ~ 2.5. 

• SMGs show rather strong correlations 

with correlation lengths larger than those 
of other high redshift sources. 

• LIRGs and ULIRGs cannot be identified 
with any of the distant populations found 
by rest-frame ultraviolet and optical sur- 
veys. 

Although these findings are answering the ba- 
sic questions about the sources that make up 
the GIB, there are still observational difficulties 
to be overcome to complete these answers. The 
SEDs of LIRGs and ULIRGs are quite variable 
and often not very well constrained in their ratio 
of far-infrared to mid-infrared or to submillime- 
ter wavelengths. The far-infrared, where most 
of the energy is radiated, requires cryogenically 
cooled telescopes. These have small diameters 
and, hence, poor angular resolution and severe 
confusion limits for blind surveys. Establish- 
ing proper SEDs for the different classes of in- 
frared galaxies detected either in mid-infrared 
(with ISOGAM at 15 ^m or MIPS at 24 /im) 
or in millimeter- submillimeter surveys is one of 
the challenges of the coming decade. Making 
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sure that no class of sources that contribute 
significantly to the CIB at any wavelength has 
been missed is an other observational challenge. 
The submillimeter galaxies not found through 
the radio-selected sources and the question of 
the warm submillimeter galaxies are also two of 
those challenges. 

Multiwavelength observations of high-2; in- 
frared galaxies give a number of new insights 
on the galaxy formation and evolution prob- 
lem. As an example, the gas masses and to- 
tal masses of SMGs are found to be very high. 
There is a first indication that the number of 
such high-mass object at redshifts between 2 
and 3 is uncomfortably largo compared to semi- 
analytical models of galaxy formation based on 
the standard hierarchical structure-formation 
frame. The evolution of the luminosity function 
is dominated by more luminous sources as red- 
shift increases. This is surprising because the 
mass function of the collapsed structure is ex- 
pected to be dominated by smaller and smaller 
objects as redshift increases. 

The populations of infrared galaxies concen- 
trated at 2 ~ 1 and at z ~ 2.5 studied so 
far reveal rather different type of sources. The 
lower redshift ones seem to be starburst phases 
of already-built massive, late- type field galax- 
ies accreting gas or gas-rich companions form- 
ing the disks. We see today a rapid decrease 
of this activity probably associated with a dry 
out of the gas reservoir in their vicinity. The 
larger redshift ones, which are also more lumi- 
nous, seem to belong to more massive complex 
systems involving major merging. These sys- 
tems could be located in the rare larger ampli- 
tude peaks of the large-scale structures leading 
to massive elliptical galaxies at the center of 
rich clusters. The redshift distribution of these 
seems quite similar to the redshift distribution 
of quasars. 

Interesting problems that are central to the 
understanding of galaxy formation and evolu- 
tion have to be solved in the next decade: 

• Determining the role of the large scale en- 
vironment (nodes, filaments and sheets of 
the large-scale structures) on star forma- 
tion; 

• Find the relative rates of accretion of gas 
and smaller galaxies in the growth of mas- 
sive objects; 



• Establish the cycle of bulge versus disk 
formation, as a function of the ratio be- 
tween stars and gas in the accreted mate- 
rial; 

• Identify the different types of starburst 
(in a disk or in the nucleus, interaction 
or merger driven); and 

• Estimate the fraction of time spent in the 
starburst phase and the duration of this 

phase 

Current observations all point in the direction 
of a possible strong effect of the large-scale envi- 
ronment and the need for models of hierarchical 
formation and evolution that include properly 
star versus gas ratio in the accreted material. 

Finally the connection between the starburst 
phenomenon and the AGN activity is an old 
question still largely unresolved. Recent obser- 
vations of infrared/submillimeter galaxies have 
reinforced the link but have not much improved 
our understanding of the physical link. 
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APPENDIX: THE PHENOMENOLOG- 
ICAL MODEL OF LAGACHE ET AL. 

(2004) 

In this paper, we make extensive use of the La- 
gache et al. (2004) phenomenological model to 
illustrate our points. This model constrains in 
a simple way the evolution of the infrared lumi- 
nosity function with redshift. It fits all the ex- 
isting source counts consistent with the redshift 
distribution, the CIB intensity, and, for the first 
time, the CIB fluctuation observations, from 
the mid-infrared to the submillimeter range. In 
this model, Lagache et al. (2004) assume that 
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infrared galaxies are mostly powered by star for- 
mation and hence they use SEDs typical of star- 
forming galaxies. Although some of the galax- 
ies will have AGN-dominated SEDs, they are a 
small enough fraction that they do not affect the 
results significantly. They therefore construct 
"normal" and starburst galaxy template SEDs: 
a single form of SED is associated with each ac- 
tivity type and luminosity. They assume that 
the luminosity function is represented by these 
two activity types and that they evolve indepen- 
dently. They search for the form of evolution 
that best reproduces the existing data. An ex- 
ample of two cosmological implications of this 
model is (a) the PAH ft^aturcs remain promi- 
nent in the redshift band 0.5 2.5 (as observa- 



tionally shown by e.g., Caputi et al. 2005), and 
{b) the infrared energy output has to be domi- 
nated by ~3'^10^^ L0 to ~3xl0^^ L© galaxies 
from redshift 0.5 to 2.5. 

The excellent agreement between the model and 
all the available observational constraints makes 
this model a likely good representation of the 
average luminosity function as a function of red- 
shift and a useful tool to discuss obsca'vations 
and models. Its rather simple assumptions such 
as the single parameter sequence of SEDs for 
starbTirst galaxies is certainly not accounting 
for some of the detailed recent observations but 
probably do not affect seriously the redshift evo- 
lution of the averaged properties which are what 
is modeled. 
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